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Q Abstract 
(N 

I The dispersion in the entire Brillouin zone and the temperature de- 

pendence (right up to the melting temperature) of the anharmonic 
frequency shift and phonon damping in a number of fee metals is 
investigated on the basis of microscopic calculations. It is found 
that the anharmonic effects depend sharply on the wave vector in 
the directions T-X, X-W, and T-L and, in contrast to bec metals, 
; the magnitude of the effects is not due to the softness of the ini- 

^ ■ tial phonon spectrum. It is shown that the relative frequency shifts 

\Q • and the phonon damping near melting do not exceed 10-20%. The 

relative role of various anharmonic processes is examined, and the 
relation between the results obtained and existing experimental data 
is discussed. 
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*^ ■ The investigation of anharmonic effects (AEs) in lattice dynamics is a classic problem of 

solid-state physics. It is important, specifically, because of the role that such effects they can 
q \ play in phenomena associated with structural phase transitions and melting in crystals (see, 
for example, Refs. 1-3). At the same time, obtaining any information about the magnitude 
and scale of AEs from experiment and theory is a difficult problem. The experimental 
study of such "basic" AEs as the frequency shift and damping of phonons is very difficult 



and leads to a large uncertainty in the results (see, for example, the data presented in 
Refs. 4 and 5 for bec and fee metals, respectively). Up to now first-principles microscopic 
calculations of AEs have been performed for one point of the Brillouin zone (N) in the bec 
phase of Zr and four points (N, P, u, G) in Mo. (Ref. 6) Detailed information about AEs 
in the entire Brillouin zone and their temperature dependence has been obtained in Refs. 4 
and 7 on the basis of pseudopotential theory for the bec phases of alkali and alkaline-earth 
metals. For these metals the most striking manifestations of AEs are due to the "soft- 
mode behavior" (the anomalous temperature dependence of the phonon frequencies) of the 
S4 branch. It is of interest to calculate AEs for the "general position", i.e., for crystals 
which do not possess soft vibrational modes. Such crystals include most metals with close- 
packed structures, for example, fee. There is virtually no information about AEs in such 
metals, not counting Ref. 5, where the damping of phonons in precious metals and Al 
was calculated. However, the calculations were performed only for individual points of the 
Brillouin zone and at room temperature (which is much less than the melting temperature 
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T m ). Moreover, the accuracy of the model used in Ref. 5 for interatomic interactions is 
unclear, which is probably the reason for the qualitative discrepancy between the theoretical 
and experimental dependences of the damping on the wave vector. In the present letter 
the AEs in the lattice dynamics of fee metals are investigated in detail on the basis of a 
systematic anharmonic perturbation theory. 

The calculations were performed for Ir, the fee phases of Ca and Sr, and the hypothetical 
fee phase of K. This choice of objects was determined by the fact that reliable and at the 
same simple microscopic models which make it possible to describe a wide range of lattice 
properties of these metals exist Refs. 4,7-9. Iridium is an example of an fee crystal with 
a "stiff interatomic interaction potential, similar to the Lennard- Jones potential, Ref. 8 
while K, Ca, and Sr are characterized by potentials with a softer "core". To demonstrate 
the characteristic features of AEs due to the specific nature of the potential with the same 
lattice geometry, we present here the results for the "limiting" cases — Ir and K. The results 
for Ca and Sr are qualitatively similar to the results for K and will be presented elsewhere. 

The parameters determining the interatomic interaction potential are presented in Ref. 
8 for Ir and in Ref. 4 for K. The calculations were performed on the basis of the anharmonic 
perturbation theory taking account of thermal expansion (quasiharmonic contributions, qh) 
and three- and four-phonon interaction processes to accuracy x 2 , where x = (m/M) 1 / 4 is 
the adiabatic smallness parameter, and m and M are the electron and ion masses. The 
exact formulas actually used in the calculations are presented in Ref. 4. For a qualitative 
discussion of the results, we present here their high-temperature asymptotic representation 
for T > Q D (Q D is the Debye temperature): 

Aau = Af k + A Ak + A 4 Ak (1) 

*ah AQ 

Af k = - lx ^—, (2) 

A 3 _Z_ \- |T,k,q,k+ q | 2 1 K 2 ~ ^I) 2 ~ ^Ak(^l + <4) /ON 

Ak - 4M3^ Ak £j | V >V» | ^2 (fa; 2 + U 2 _ ^2 _ 4^2 > W 

rAk = E K k ^ k+q 2 [<w - "i - <"2) + + - u, 2 )} (5) 
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Here A, /i, and v are the indices of the phonon branches, k and q are quasimomenta, cj Ak , 
Aaic, and Taic are, respectively, the initial phonon frequency and its shift and the phonon 
damping, V and W are the amplitudes of three- and four-phonon processes (see Ref. 4), 
7Ak = — duxk/dQ, is Gruneisen parameter. AQ is the change in the volume Q per atom due 
to thermal expansion, and 1 = (/i, q) and 2 = (v, k + q). 

The basic computational results are presented in Figs. 1-5. It is evident from Figs. 1 
and 2 that the high-temperature asymptotic representations (3)-(5) obtain very early (for 
T > 0d/3) Comparing the analogous results for bec metals 4 it can be concluded that this 
property is probably quite general. It follows from Figs. 1-4 that the strongest AEs in Ir 
occur for longitudinal phonons near the point X. Specifically, the decrease of the frequency 
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and the increase of the damping with temperature are greatest for these phonons The 
situation differs sharply in this respect from bcc alkali and alkaline-earth metals, Refs. 4,7 
where the maximum damping obtains for "soft" phonons, for which the opposite behavior 
of the frequency with increasing T is characteristic — hardening instead of softening. It is 
evident from Fig. 5 that this difference is due to the fact that in Ir three-phonon processes 
dominate over four-phonon processes. 

To understand the role of the characteristic features of the interatomic interactions for 
the same lattice geometry, the computational results tor a hypothetical fee phase of K 
are presented in Figs. 3 and 4 The explicit domination of three-phonon over four phonon 
processes, which causes the frequency shift to be negative over the entire Brillouin zone is 
also seen in this case and is apparently characteristic for all metals with fee structure At 
the same time, the fact that the neighborhood of the point X is distinguished probably 
occurs only for crystals with a "hard" interatomic interaction potential of the type in Ir. 

It is well known Ref. 8 that the "hardness" of the potential in metals increases with 
the effective valence Z (Z= 1 for K and Z=4.5 for Ir). In any case Ir is a striking example 
of a crystal where AEs are maximum for high-frequency vibrations, in contrast to crystals 
with soft modes where the AEs are maximum for low-frequency phonons. The latter can 
be easily explained by the high powers of the frequency in the denominators of the Eqs. 
(3)- (5), but the example of Ir shows that in the general case the question cannot be solved 
on the basis of such simple considerations. 

According to Fig. 4, an important feature of phonon damping is its nonmonotonic 
dependence on the wave vector in the V — L and V — X directions. According to the 
experimental data presented in Ref. 5, such nonmonotonic behavior is observed in Al, Cu, 
Ag, and Au. It can therefore be supposed that this behavior is typical for all metals with 
the fee structure. 

In closing we note that, as follows from Figs. 3-5, the AEs have a strong dependence 
over the Brillouin zone. For this reason, approximations which take them into account 
"on average" (approximations of the type used in the recently published Ref. 10), are 
dangerous. Finally, from the fact, found in the present work, that three-phonon processes 
predominate over four-phonon processes it follows that the well-known approaches such as 
the self-consistent phonon approach are inapplicable to metals with a close-packed structure 
Ref. 11. 
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Figure captions 

1. Temperature dependence of the frequency shift in Ir at the symmetric points of the 
Brillouin zone; T m — melting temperature, Qd — Debye temperature. The indices 
1 and 2 denote longitudinal and transverse phonons; uj p \ — ion plasma frequency 
{ujr = 871 K, see Ref. 8). 

2. Temperature dependence of phonon damping in Ir. All notations are the same as in 
Fig 1. 

3. Dispersion of the relative frequency shift 5 = A/u at T = T m . Solid line - Ir, dotted 
line - fee K. The numbers 1, 2 indicate longitudinal and transverse branches. 

4. Dispersion of the relative phonon damping r\ = T/u at T = T m . Solid line - Ir, dotted 
line - fee K. The numbers 1, 2 indicate longitudinal and transverse branches. 

5. Variation of the anharmonic contributions A qh (a), A 3 (b), and A 4 (c) in the T—X—W 
directions in Ir at T = T m . 
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